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Abstract 

A search is presented for photonic signatures motivated by generalized models of gauge- 
mediated supersymmetry breaking. This search makes use of 20.3 fb -1 of proton-proton col¬ 
lision data at = BTeV recorded by the ATLAS detector at the LHC, and explores models 
dominated by both strong and electroweak production of supersymmetric partner states. Four 
experimental signatures incorporating an isolated photon and significant missing transverse 
momentum are explored. These signatures include events with an additional photon, lepton, 
&-quark jet, or jet activity not associated with any specific underlying quark flavor. No signific¬ 
ant excess of events is observed above the Standard Model prediction and model-dependent 
95% confidence-level exclusion limits are set. 
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1 Introduction 

This paper reports on a search for four classes of events containing energetic isolated photons and large 
missing transverse momentum (with magnitude denoted £™ ss ) in 20.3 fb~ 1 of proton-proton (pp) collision 
data at Vs = 8 TeV recorded with the ATLAS detector at the Large Hadron Collider (LHC) in 2012. 
For the first of the four classes, two isolated energetic photons are required (“diphoton” events), while 
for the remaining classes only a single isolated photon is required. For the second and third classes, 
the isolated photon is required to appear in combination with a “b- jet” identified as having arisen from 
the production of a bottom ( b ) quark (“photon+Z?” events) or an isolated electron or muon (“photon+r’ 
events), respectively. For the fourth class of events the isolated photon is required to appear in combination 
with multiple jets selected without regard to the flavor of the underlying parton (“photon+/’ events). 

The results are interpreted in the context of a broad range of general models of gauge-mediated su¬ 
persymmetry breaking (GGM) [1-3] that include the production of supersymmetric partners of strongly 
coupled Standard Model (SM) particles as well as SM partners possessing only electroweak charge. In 
all models of GGM, the lightest supersymmetric particle (LSP) is the gravitino G (the partner of the 
hypothetical quantum of the gravitational field), with a mass significantly less than 1 GeV. In the GGM 
models considered here, the decay of the supersymmetric states produced in LHC collisions would pro¬ 
ceed through the next-to-lightest supersymmetric particle (NLSP), which would then decay to the G LSP 
and one or more SM particles, with a high probability of decay into y + G. In this study, several different 
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possibilities for the nature of the NLSP are considered, providing separate motivation for the four different 
and complementary experimental signatures that are explored. In all models considered, all supersymmet¬ 
ric states with the exception of the G are short lived, leading to prompt production of SM particles that are 
observed in the ATLAS detector. 

The results based on the diphoton and photon+b signatures extend and supplant studies (Refs. [4] 
and [5], respectively) that made use of 4.8 fb~ 1 of pp collision data at \[s = 1 TeV; the analyses based on 
the photon +j and photon+7 signatures are new and have only been performed with the 8 TeV data. Making 
use of 19.7 ftr'of pp collision data at yfs = 8 TeV, a search [6] for events similar in nature to those of the 
diphoton and photon+ j signatures mentioned above has performed by the CMS Collaboration, and used 
to set limits on the masses of strongly coupled supersymmetric particles in several GGM scenarios. 

2 Gauge-mediated supersymmetry phenomenology 

Supersymmetry (SUSY) [7-15] introduces a symmetry between fermions and bosons, resulting in a SUSY 
partner (sparticle) with identical quantum numbers except a difference by half a unit of spin for each SM 
particle. As none of these sparticles have been observed, SUSY must be a broken symmetry if realized 
in nature. Assuming R-parity conservation [16-20], sparticles are produced in pairs. These would then 
decay through cascades involving other sparticles until the stable, weakly interacting LSP is produced, 
leading to a final state with significant E ™ 1SS . 

Experimental signatures of gauge-mediated SUSY breaking models [21-26] are largely determined 
by the nature of the NLSP. For GGM, the NLSP is often formed from an admixture of any of the SUSY 
partners of the electroweak gauge and Higgs boson states. In this study, three cases are assumed for the 
composition of the NLSP. For the first case, the NLSP is assumed to be purely binolike [the SUSY partner 
of the SM U(l) gauge boson]. For the second case, the NLSP is assumed to be an admixture of bino 
and neutral higgsino states. For the final case, the NLSP is assumed to be a degenerate triplet of wino 
states [the SUSY partners of the SM SU(2) gauge bosons]. In this paper, the neutral NLSP is denoted^ 
irrespective of its composition. For the case that the NLSP is a degenerate triplet, the charged NLSP states 
are denoted^*. The properties of the GGM models used to represent these possibilities are discussed 
below and summarized in Table 1 . 

For the case that the NLSP is a bino, the final decay in each of the two cascades in a GGM event 
would be predominantly x® —» y + G, leading to final states with yy + /il [ n ' ss . For the case that the NLSP 
is a mixture of the bino and higgsino, both the possibilities that the higgsino mass parameter p is less 
than or greater than zero are explored. For the p < 0 possibility, the final decay in the cascade would 
include a significant contribution from x® -* h + G with the subsequent decay h —> bb, leading to final 
states with a photon, multiple Z?-jets, and Z7™ ss . The latter (p > 0) possibility can produce scenarios for 
which the final decay in the cascade can be relatively evenly split between x® —> y + G and —> Z + G, 
leading to final states with a photon, multiple jets (including two from the hadronic decay of the Z boson) 
that most often do not arise from ^-quarks, and /7. l [ n ' ss . For the case that the NLSP is a degenerate set of 
three wino states, the final step in the cascade includes charged as well as neutral wino decays. Charged 
wino decays tend to produce isolated leptons, while neutral wino decays produce photons with a wino- 
to-photon branching fraction that is no less than sin 2 8\y for any value of the wino mass. Overall, these 
two wino-NLSP contributions lead to a significant number of events with an isolated photon accompanied 
by an isolated lepton. Of the five GGM models considered here, two (the “gluino-bino” and “wino-bino” 
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Figure 1: Typical production and decay-chain processes for the gluino-production (left) and electroweak- 
production (right) instances of the GGM model for which the NLSP is a binolike neutralino, referred to in 
the text as the gluino-bino and wino-bino models, respectively. 


models, where the gluino is the SUSY partner of the gluon) incorporate a purely binolike NLSP, two 
(the “higgsino-bino” models) incorporate a NLSP that is a higgsino-bino admixture, and one (the “wino- 
NLSP” model) incorporates a winolike set of NLSPs; in all cases the mass of the NLSP state is considered 
to be a free parameter of the model. 

The two GGM models incorporating a binolike NLSP are the focus of the diphoton analysis. For these 
models, one other set of SUSY partner states is taken to be potentially accessible in 8 TeV pp collisions, 
while all other SUSY masses are decoupled (set to inaccessibly large values). For both of these binolike 
NLSP cases, production proceeds solely through this set of SUSY partners, with the NLSP appearing in 
the subsequent decays of the produced SUSY partner states. For the gluino-bino model, the set of partners 
is composed of a degenerate octet of gluinos. For the wino-bino model, the set of partners is composed of 
a degenerate triplet of wino states x* 2 >X\> and is dominated by the production ofx\ x\ and;^ xf- For both 
of these models, the masses of these produced states are considered to be free parameters along with that 
of the chosen^ state, the latter of which is constrained to be less than those of the produced states. This 
results in a SUSY production process that proceeds through the creation of pairs of the higher-mass states, 
which subsequently decay through short cascades to the NLSPjf^ states. Other SM objects (jets, leptons, 
photons) may be produced in these cascades. The branching fraction to y + G is 100% for m~a —> 0 
and approaches cos 2 dw for m~n » m z , with the remainder of the sample decaying to Z + G. For all 
X° { masses, then, the branching fraction is dominated by the photonic decay, leading to the diphoton-plus- 
Z?™ ss signature. For these models with a binolike NLSP, typical production and decay channels for strong 
(gluino) and electroweak (wino) production are exhibited in Fig. 1. 

The higgsino-bino GGM models incorporate a NLSP composed of a higgsino-bino admixture, as well 
as a degenerate octet of gluinos identical in nature to those of the gluino-bino model. For the first of these 
models, which is the focus of the photon+fi analysis, the higgsino mass parameter p is required to be 
negative, and the composition of the NLSP is set by adjusting p and the GGM U(l) mass parameter M\ 
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Table 1: Summary of the five GGM models considered in this study. For the two higgsino-bino models, 
the functions f±(M\,n ) are chosen to establish NLSP decay properties commensurate with the target 
experimental signature, as described in the text. 


GGM Model 

Experimental 

Signature 

Produced 

State(s) 

Composition 
of NLSP 

Free 

Parameters 

Gluino-bino 

diphoton 

gluino 

bino 

M Zj, M x o 

Wino-bino 

diphoton 

wino 

bino 


Higgsino-bino (/r < 0) 

photon+6 

gluino, higgsino 

higgsino/bino 

MgJ.(M u n) 

Higgsino-bino (ju > 0) 

photon+ /' 

gluino, higgsino 

higgsino/bino 

Mg,f + (M u /u) 

Wino NLSP 

photon+f 

wino 

wino 

M w 



Figure 2: Typical production and decay-chain processes for the gluino-production instance of the GGM 
model for which the NLSP is a higgsino-bino neutralino admixture, referred to in the text as the higgsino- 
bino model. For the model with fi < 0 (left), the final step of the cascade (the x® decay) would have a 
probability of order 50% of producing a Fliggs boson rather than a photon or Z boson; for the model with 
H > 0 (right), the x® decay would have a probability of order 50% of producing a Z boson rather than 
a photon. For both of these models, production can also proceed through gaugino and neutralino states, 
which can dominate the production cross section for high values of gluino mass. 
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Figure 3: Typical production and decay-chain processes for the wino-NLSP model. In this model, the x® 
is a pure W° state, while the x\ are the two charged wino states. 

so that a constant ratio of the branching fraction of x~ —> h + G to that of jj>° —> y + G is maintained at 
approximately 1.7:1 over the full range of NLSP masses. The photon+fi analysis was found to provide 
the greatest advantage relative to the diphoton analysis for this ratio of branching fractions. In the limit 
that m -o » mz, the NLSP branching fractions to h + G,y + G, and Z + G approach 56%, 33%, and 11%, 
respectively. The GGM SU(3) mass parameter M 3 bears a direct relation to the gluino mass, and is taken 
to be a free parameter in this n < 0 higgsino-bino model, with all squark states decoupled. The GGM 
SU(2) mass parameter M 2 is set to a value of 2.5 TeV. Four other electroweak gaugino states typically lie 
within 25 GeV of the x® NLSP: the two lightest charginos x±, and two additional neutralinos Xn and 
The pair production of gluinos or any of these four additional gaugino states leads to decays to the x® via 
cascades involving SM particles. 

For the second of the higgsino-bino models, which is the focus of the photon +j analysis, the p para¬ 
meter is chosen to be positive, which suppresses the h + G decay mode of the higgsino. As in the models 
described above, the NLSP mass is taken to be a free parameter. The M\ and // parameters are adjusted so 
that the branching fractions of thejf-j to y + G.Z + G and h + G are maintained close to 50%, 49% and 1% 
for most values of thejj^ and gluino masses. In this model, the production of gluino pairs can be followed 
by decays to both a single photon and a hadronically decaying Z boson, producing events with a single 
isolated high-energy photon accompanied by two jets. In the case that the gluino mass is substantially 
larger than the X® mass, additional jets can be produced in the cascade. Three additional electroweak 
gaugino states lie close in mass to the x®, allowing for the possibility of SUSY production through pairs 
of these states. Such events tend to produce fewer jets than those that proceed through gluino production, 
but in certain regions of the model space can provide a significant contribution to data samples selected 
to isolate the photon-plus-jets signature. As in the /r < 0 higgsino-bino model, the value of M 3 , which is 
directly related to the gluino mass, is taken to be a free parameter, M 2 is set to a value of 2.5 TeV, and all 
squark states are decoupled. Typical production and decay-chain processes for the two models for which 
the NLSP is a higgsino-bino admixture are shown in Fig. 2. 

Finally, the wino-NLSP model, which is the focus of the photon+f analysis, incorporates a set of three 
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degenerate winolike NLSPs. This set includes the neutral IT 0 , which as the lightest neutral gaugino is 
also referred to as the X®, as well as the two charged wino states, which form the x\ states. Production 
proceeds through the direct production of pairs of NLSP states; such events usually contain at least one 
W° NLSP Although the VF° couples preferentially to the Z boson relative to the photon, the W° decays 
into a photon+gravitino final state with unit branching fraction for wino mass below that of the Z boson. 
The W° branching fraction to photon+gravitino approaches sin 2 %- for wino masses far above that of the 
Z boson. Leptons can be produced either through the decays of charged wino states, or through the decays 
of Z bosons that arise from W° decay, leading to a significant probability that the overall final state would 
contain both a photon and a lepton. In this model, a common wino mass scale is taken as a free parameter, 
with all other GGM mass parameters set to a value of 2.5 TeV, except the squark masses, which are set to 
infinity. A production and decay diagram typical for this model is shown in Fig. 3. 

For all five models considered here, the mass of the gravitino is chosen so that the NLSP decay length is 
never greater than 1 mm. This ensures that all particles arising from the decay of the NLSP are prompt, and 
in particular that the relationship between the point and direction of impact of photons from NLSP decay 
upon the face of the detector is consistent with that of a prompt photon (a separate analysis [27] searches 
for GGM models with a longer-lived binolike NLSP, leading to signatures with nonprompt photons). In 
addition, the ratio tan/3 of the two SUSY Higgs-doublet vacuum-expectation values is set to a value of 
1.5; for all five models, the phenomenology relevant to this search is only weakly dependent on the value 
of tan /3. 

3 Samples of simulated processes 

For the GGM models under study, the SUSY mass spectra and branching ratios are calculated using 
SUSPECT 2.41 [28] and SDECAY 1.3b [29], respectively, inside the package SUSY-HIT 1.3 [30]. The 
Monte Carlo (MC) SUSY signal samples are produced using HERWIG++ 2.5.2 [31] with CTEQ6L1 parton 
distribution functions (PDFs) [32]. Signal cross sections are calculated to next-to-leading order (NLO) 
in the strong coupling constant, including, for the case of strong production, the resummation of soft 
gluon emission at next-to-leading-logarithmic accuracy (NLO+NLL) [33-37], The nominal cross section 
and its uncertainty are taken from an envelope of cross-section predictions using different PDF sets and 
factorization and renormalization scales [38]. At fixed center-of-mass energy, SUSY production cross 
sections decrease rapidly with increasing SUSY partner mass. At yfs = 8 TeV, the gluino-production 
cross section is approximately 24 fb for a gluino mass of 1000 GeV and falls to below 1 fb for a gluino 
mass of 1400 GeV. The wino-production cross section is approximately 15 fb for a wino mass of 500 GeV, 
and falls to approximately 1 fb for a wino mass of 750 GeV. 

While most of the backgrounds to the GGM models under examination are estimated through the use 
of control samples selected from data, as described below, the extrapolation from control regions (CRs) 
to signal regions (SRs) depends on simulated samples, as do the optimization studies. The simulation 
of W and Z boson production, including events with up to five accompanying partons, is calculated by 
two different generators. The ALPGEN 2.14 [39] Monte Carlo generator is interfaced to HERWIG 6.520 
for showering and fragmentation and to JIMMY [40] for simulation of the underlying event. Parton dis¬ 
tributions are provided by the CTEQ6L1 functions. Similar samples are produced with the SHERPA 1.4.1 
generator [41] with CT10 [42] PDFs, for up to four accompanying partons. 

Wy production is also simulated via ALPGEN interfaced to HERWIG and JIMMY, but makes use of the 
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CT10 PDFs. Other Wy samples are generated, as is the Zy process, by using SHERPA with the CT10 
PDFs. The tty process is simulated at leading order (LO) using MADGRAPH 5.1.5.11 [43] and CTEQ6L1, 
interfaced to the PYTHIA 6.427 parton shower generator [44], The tt process is simulated not only with 
the the POWHEG generator interfaced to PYTHIA and the CTEQ6L1 PDFs, but also with the MCONLO 4.06 
generator [45,46] and the CT1Q PDFs, including full NLO QCD corrections. This contribution is rescaled 
to match the tt cross section at NNLO with NNLL soft gluon terms, as calculated with top++2.0 [47-52]. 
The ty and ty processes are simulated with the WHIZARD 2.1.1 [53,54] generator, with four-flavor/five- 
flavor matching provided using HOPPET [55]. Additional photon radiation is added with PHOTOS [56], with 
parton showering and fragmentation again simulated with PYTHIA. Other t and t samples are generated 
with POWHEG. 

The y+jet(s) process is simulated in a similar manner to the W ± or Z samples using ALPGEN interfaced 
to HERWIG and JIMMY and the CTEQ6L1 PDFs. A generator-level requirement of 35 GeV is applied to the 
photon transverse momentum Pj, and the sample is generated in exclusive bins of p y T to produce a more 
statistically significant sample at higher values of p^. Additional y+jet(s) samples are used, simulated with 
SHERPA and the CT10 PDFs. The prompt diphoton sample is generated with PYTHIA 6.423, which includes 
the subprocesses gg —» yy and qq —* yy, with the requirement that there be at least two prompt photons 
with generated transverse momentum greater than 20 GeV. Parton densities are modeled according to the 
MRST 2007 LO* [57] functions. 

The background from Z(—> vv) + yy production is simulated using the SHERPA MC generator, normal¬ 
ized to a cross section calculated at LO using MADGRAPH 5 and the CTEQ6L1 PDF, and then corrected by a 
V-factor of 2.0 ±1.0 [58]. The background from W (—> iv) + yy production is simulated using the ALPGEN 
MC generator, although the overall normalization is set via a study making use of data events containing 
two photons and a charged lepton (to be discussed below). Diboson production, for the case that each 
boson is a IT or Z, is simulated with POWHEG. 

All MC samples are processed with the GEANT4-based simulation [59] of the ATLAS detector [60], or, 
where appropriate, a simulation of the ATLAS detector based on parametrized shower shapes in the calor¬ 
imeter, and GEANT4 elsewhere. Corrections are applied to the simulated samples to account for differences 
between data and simulation for the lepton and photon trigger, identification, and reconstruction efficien¬ 
cies, as well as for the efficiency and misidentification rate of the algorithm used to identify jets containing 
/;-hadrons (^-tagging). The variation of the number of pp interactions per bunch crossing (“pileup”) as 
a function of the instantaneous luminosity is taken into account by overlaying simulated minimum-bias 
events according to the observed distribution of the number of pileup interactions in data, with an average 
of 21 interactions per event. 

4 ATLAS detector 

The ATLAS experiment makes use of a multipurpose detector [61] with a forward-backward symmetric 
cylindrical geometry and nearly 4/r solid angle coverage. 1 Closest to the beam line are solid-state tracking 


'ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the 
detector and the z axis along the beam pipe. The x axis points from the IP to the center of the LHC ring, and the y axis points 
upwards. Cylindrical coordinates (r, <j>) are used in the transverse plane, (f> being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as 77 - - ln[tan(#/2)]. Angular distance is measured in units of 
A R= V(At; ) 2 + (A <t>f. 
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devices comprising layers of silicon-based pixel and strip detectors covering |t/| < 2.5 and straw-tube 
detectors covering |//| < 2.0, located inside a thin superconducting solenoid that provides a 2T magnetic 
field. Outside the solenoid, fine-grained lead/liquid-argon electromagnetic (EM) calorimeters provide 
coverage over |//| <3.2 for the measurement of the energy and position of electrons and photons. A 
presampler, covering \r/\ < 1.8, is used to correct for energy lost upstream of the EM calorimeter. An 
iron/scintillator-tile hadronic calorimeter covers the region \r/\ < 1.7, while a copper/liquid-argon medium 
is used for hadronic calorimeters in the end cap region 1.5 < |//| < 3.2. In the forward region 3.2 < 
|//| < 4.9 liquid-argon calorimeters with copper and tungsten absorbers measure the electromagnetic and 
hadronic energy. A muon spectrometer consisting of three superconducting toroidal magnet systems, each 
comprising eight toroidal coils, tracking chambers, and detectors for triggering, surrounds the calorimeter 
system. The muon system reconstructs penetrating tracks over a range |//| < 2.7 and provides input to the 
trigger system over a range |//| < 2.4. A three-level trigger system is used to select events. The first-level 
trigger is implemented in hardware and uses a subset of the detector information to reduce the accepted 
rate to at most 75 kHz. This is followed by two software-based trigger levels that together reduce the 
accepted event rate to 400 Hz on average depending on the data-taking conditions during 2012. 

5 Reconstruction of candidates and observables 

Primary vertices are formed from sets of two or more tracks, each with transverse momentum p^ ack > 
400 MeV, that are mutually consistent with having originated at the same three-dimensional space point 
within the luminous region of the colliding proton beams. When more than one such primary vertex is 
found, the vertex with the largest scalar sum of the squared transverse momenta of the associated tracks 
is chosen. To further ensure the event resulted from a beam collision, the primary vertex of the event is 
required to have at least five associated tracks. 

Electron candidates are reconstructed from EM calorimeter energy clusters consistent with having 
arisen from the impact of an electromagnetic particle (electron or photon) upon the face of the calorimeter. 
For the object to be considered an electron, it is required to match a track identified by a reconstruction 
algorithm optimized for recognizing charged particles with a high probability of bremsstrahlung. In ad¬ 
dition, the matched track is required to include information from at least seven layers of the solid-state 
tracking system; a track within the acceptance of the tracking system typically traverses eleven layers of 
the solid-state tracking system. The energy of the electron candidate is determined from the EM cluster, 
while its pseudorapidity is determined from the associated reconstructed track. Further details of the re¬ 
construction of electrons can be found in Refs. [62] and [63]. Electron candidates used by these analyses 
are further required to have p \ > 20GeV and \rj\ < 2.47. For the photon+( analysis, signal electrons are 
not allowed to be within the transition region 1.37 < \r/\ < 1.52 between the barrel and end cap calori¬ 
meters. A track-based isolation requirement is imposed, with the scalar sum of the transverse momenta of 
tracks within a cone of size A R = 0.3 required to be less than 16% of the electron p-\ . Finally, the electron 
track is required to be consistent with coming from the primary vertex in the r-z plane. 

Electromagnetic clusters are classified as photon candidates provided that they either have no matched 
track or have one or more matched tracks consistent with coming from a photon conversion vertex. 
Based on the characteristics of the longitudinal and transverse shower development in the EM calori¬ 
meter, photons are classified as “loose” or “tight.” Further details of the reconstruction of photons can be 
found in Ref. [64]. In the case that an EM calorimeter deposition is identified as both a photon and an 


electron, the photon candidate is discarded and the electron candidate retained. Photon candidates used by 
these analyses are required to be within \p\ <231, and to be outside the transition region 1.37 < |//| < 1.52. 
Finally, an isolation requirement is imposed. After correcting for contributions from pileup and the de¬ 
position ascribed to the photon itself, loose and tight isolation criteria are defined, with the tight criterion 
requiring less than 4 GeV of transverse “isolation energy” in a cone of size A R = 0.4 surrounding the 
energy deposition in the calorimeter associated with the photon. For the loose isolation criterion, no more 
than 5 GeV of isolation energy is allowed within a cone of size A R = 0.2. The tight criterion is used 
for the diphoton analysis, while the loose criterion is used for the remaining three signatures (photon+ 6 , 
photon+ j, photon+f). 

Muon candidates make use of reconstructed tracks from the tracking system as well as information 
from the muon system [65]. Muons are required to be either “combined,” for which the muon is recon¬ 
structed independently in both the muon spectrometer and the tracking system and then combined, or 
“segment-tagged,” for which the muon spectrometer is used to tag tracks as muons, without requiring a 
fully reconstructed candidate in the muon spectrometer. Signal muons are required to have pj > 20 GeV 
and |;/| < 2.4. Track-based as well as calorimeter-based isolation requirements are imposed, with the 
scalar sum of the transverse momenta of tracks within a cone of size A R - 0.3 required to be less than 
12 % of the muon p-\, and the energy in the calorimeter projected in the transverse plane within a cone of 
size A R = 0.3, corrected for pileup, also required to be less than 12% of the muon p-\. Finally, the muon 
track is required to be consistent with coming from the primary vertex in both the r-z and r-(p planes. 

Jets are reconstructed from three-dimensional calorimeter energy clusters using the anti-A algorithm [ 66 ] 
with a radius parameter R = 0.4. Jets arising from detector noise, cosmic rays or other noncollision sources 
are rejected, as described in Ref. [67]. Each cluster is classified, prior to the jet reconstruction, as com¬ 
ing from an electromagnetic or hadronic shower on the basis of its shape [ 68 ]. Each cluster energy is 
then corrected by weighting electromagnetic and hadronic energy deposits with correction factors derived 
from Monte Carlo simulation. A correction is applied to subtract the expected contamination from pileup, 
calculated as the product of the jet area in 77—0 space and the average energy density of the event [69]. 

A further calibration, relating the response of the calorimeter to in situ jet-energy measurements [69] is 
then applied. Once calibrated, jets are required to have pj > 20GeV and |//| < 2.8. Jets containing 6 - 
hadrons are identified using the MVlc 6 -tagging algorithm [70], This neural network algorithm combines 
the information from various algorithms based on track impact-parameter significance or explicit recon¬ 
struction of b- and c-hadron decay vertices. The analyses presented in this paper use an operating point 
corresponding to 70% efficiency for jets originating from the fragmentation of a 6 -quark in simulated tt 
events, selecting approximately 0.7% of light-quark and gluon-induced jets and 20% of c-quark-induced 
jets. 

In the case that two reconstructed objects are in close enough proximity to one another to raise a 
concern that they are a single detector object reconstructed as more than one particle or jet candidate, an 
overlap-removal procedure is followed. To reduce the rate of electrons misidentified as photons, if the 
angular distance A R between a reconstructed electron and photon is less than 0 . 01 , the object is classified 
as an electron. 

To avoid ambiguity that arises when an electron or photon is also reconstructed as a jet, if a jet and 
an electron or photon are reconstructed within an angular distance A R = 0.2 of one another, the electron 
or photon is retained and the jet is discarded; if 0.2 < A R < 0.4 then the jet is retained and the electron 
or photon is discarded. Finally, in order to suppress the reconstruction of muons arising from showers 
induced by jets, if a jet and a muon are found with A R < 0.4 the jet is retained and the muon is discarded. 
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The vector momentum imbalance in the transverse plane is obtained from the negative vector sum 
of the reconstructed and calibrated physics objects and is referred to as missing transverse momentum 
E-'|i lss [71]. Calorimeter energy deposits arc associated with a reconstructed and identified high-px object 
in a specific order: electrons with pj > 10 GeV, photons with pj > 10 GeV, and jets with pj > 20 GeV. 
Deposits not associated with any such objects are also taken into account in the E! [ riISS determination, as 
are muons with pj > 10 GeV. 

The transverse mass Mj of a system of two massless particles with four-vectors p\ and pi is given by 

Mj = -^2p r Y\pjp(l — cos A(pip), 

where Acpi^ is the angular separation between the two vectors projected into the transverse plane. The 

y^E miss 

analyses presented here make use of the transverse mass of both the photon-E™ ss ( M. y T ) and lepton- 

£ j^miss 

E™ ss ( M. y 1 ) systems, where the lepton is taken to be massless in the transverse-mass determination. 

Several additional observables are defined to help in the discrimination of SM backgrounds from 
potential GGM signals. The total visible transverse energy Hj is calculated as the scalar sum of the 
transverse momenta of the selected photons and any additional leptons and jets in the event; a similar 
observable based only on the momenta of jets in the events is referred to as H^ ts . The “effective mass” 
772 e ff is defined as the scalar sum of Hj and £™ ss . The photon-£™ lss separation A0(y, E™ 88 ) is defined 
as the azimuthal angle between the missing transverse momentum vector and the selected photon. In 
the case of the diphoton analysis, Ar/> m ; n (y, £™ ss ) is defined to be the minimum value of AA(y, £™ 88 ) of 
the two selected photons. The minimum jct-E™ 88 separation A0 m j n (jet, £™ 188 ) is defined as the minimum 
azimuthal angle between the missing transverse momentum vector and the leading (highest-px) jets in the 
event. The number of leading jets used differs depending on the signature under study and is shown in 
Tables 2 and 3. For the diphoton analysis, leading jets are required to have px > 75 GeV, and if no such jet 
is found, no requirement is placed on the observable. The quantity A</> m j n (jet, y) is defined as the minimum 
separation between the selected photon and each of the two leading jets in the event. The quantity A R{£, y) 
is defined as the distance in rj-<p space between the leading photon and lepton. Finally, the quantity R 4 is 
defined as the scalar sum of the transverse momentum of the four highest-px jets in the event divided by 
the sum of the transverse momentum of all jets in the event. 

6 Event selection 

The data sample is selected by a trigger requiring the presence of one loose photon with energy projected 
into the plane transverse to the beam pipe (Ex) of greater than 120 GeV for the photon+/?, photon-i -j 
and photon+7 analyses, or two loose photons with Ej > 40 GeV for the diphoton analysis. Events are 
removed from the data sample if they contain jets likely to be produced by beam backgrounds, cosmic 
rays or detector noise, as described in Ref. [67], After applying data-quality requirements related to the 
beam and detector conditions, the total available integrated luminosity is 20.3 fb 1 . The uncertainty on 
the integrated luminosity is ±2.8%, estimated via the methodology of Ref. [72]. 

For the diphoton analysis, geared towards the exploration of the gluino-bino and wino-bino GGM 
models incorporating a purely binolikcy-p two separate SR selection strategies were developed: a “SR 77 ” 
selection geared towards the production of higher-mass strongly coupled SUSY states (gluinos and squarks) 
and a “SR^” selection geared towards the production of lower-mass weakly coupled SUSY states (wi¬ 


ld 



nos). For each of these approaches, two SRs are defined: the first (SR^7 L , SR^J L ) optimized for the case 
of a lower-mass x® and the second (SR^7 H , SR^_ H ) for a higher-mass;^. 

For the photon+fi analysis, geared towards the higgsino-bino GGM model with a negative value of the 
p parameter, two SRs (SRj /b , SRj' b ) are defined. The SRs are again distinguished by their optimization for 
low and high y* 1 mass, respectively. In particular, the SRj lb selection is designed to have a high acceptance 
for events that arise through the production of pairs of weakly coupled SUSY partners, which can have 
a significant cross section for the low-y^-mass reaches of the higgsino-bino GGM model explored here. 
For the photon-i- j analysis, geared towards the higgsino-bino GGM model with a positive value of the p 
parameter, a further set of two SRs are defined (SRj /J , SR^). These two SRs are once again distinguished 
by their optimization for low and high y* 1 mass, respectively. 

A final “SR^J ” signal region was developed to search for phot on+7 events arising from the GGM 

model with a winolike set of NLSPs. This SR is divided into two subsets—SR^ f and SRj /—depending 
on the flavor of the leading lepton (electron or muon). 

All four diphoton SRs require two tight, isolated photons with Ej > 75 GeV, while the SR^ b and SRj /b 
signal regions require a single tight, isolated photon with E\ > 125 GeV and Ej > 150 GeV, respectively, 
and the SRj J and SRj'j 1 signal regions require a single tight, isolated photon with Ej > 125 GeV and 
Ej > 300 GeV, respectively. The SR^ signal region requires a single tight, isolated photon with Ej > 
125 GeV. Along with £™ 1SS , leptonic, and (Z?-)jet activity, requirements are made on a number of additional 
observables, with values chosen to optimize the sensitivity to the GGM signal of interest for each SR. To 
ensure that the 72 l [ niss observable is accurately measured, minimum requirements on A(/> mm (y,is™ ss ) and 
A^min(j e fi £™ ss ) are considered for each SR. For the SR^ signal region of the photon+ j analysis, rejecting 
events with jets misidentified as photons by placing a requirement on Ar/; mm (jct, y) is found to improve the 
sensitivity of the analysis. 

To exploit the high-energy scale associated with SUSY production at masses close to the expected 
limit of sensitivity of the various SRs, several SRs include minimum requirements on one of the two 
total-transverse-energy observables Hj or m e As an illustration. Fig. 4 (left) shows the m e fi distribution 
of selected diphoton events as well as that expected from several SM sources and from characteristic 
strong-production points of the binolike NLSP GGM model. For electrons from W boson decay that 

£miss 

are misreconstructed as photons, the transverse mass Mj T of the photon-E) [ riISS system tends to be less 
than that of the W boson; because of this, the photon+/; analysis is found to benefit from a minimum 

y ^miss k 

requirement on Mj T . The SR, analysis also benefits from a requirement that the invariant mass 
of the system formed by the two most energetic fi-jets be close to the Higgs boson mass. A minimum 

£ ^miss 

requirement on the transverse mass M T T of the lepton-E™ ss system is similarly found to be effective 
in rejecting backgrounds from W boson and semileptonic tt decay for the photon+7 analysis. A further 
requirement that the electron-photon system invariant mass not be close to the Z boson mass helps to reject 
Z boson backgrounds to the photonanalysis. A requirement that SRj lb signal events have no identified 
charged leptons helps to reduce the background from semileptonic tt events, while a requirement that SR^j 
signal events have Rj < 0.85 helps reduce the background from SM events, which tend to have fewer and 
softer jets than do signal events; as an illustration, see Fig. 4 (right). Finally, a requirement that the total 
transverse energy from jets with p-y > 40 GeV (Hj ts ) be less than 100 GeV helps reduce the backgrounds 
to SR^ due to top quark production. 
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Figure 4: (Left) Distribution of m e ff, the sum of the total visible transverse energy and £™ ss , for selected 
diphoton events, after requiring A0 m i n (jet, fi™ ss ) > 0.5 but before application of a requirement on E™ ss 
and Ad> m ; n (y, E™ 1SS ). Also shown are the expected contributions of SM processes, estimated as described 
in Sec. 7, as well as the expected m e s distributions for the o) = (1300,150) GeV and (m?,, m D o) = 

y X\ J X\ 

(1300,1050) GeV gluino-bino GGM models. (Right) Distribution of R 4 , the scalar sum of the transverse 
momentum of the four highest-/; 7 - jets in the event divided by the sum of the transverse momentum of 
all jets in the event, for the sample surviving all SR” selection criteria except the R 4 { requirement itself. 
Also shown are the expected contributions of SM processes, estimated as described in Sec. 7, as well as 
the signal expectation for the two points in the M 3-/1 parameter space characteristic of the // > 0 GGM 
model relevant to the photon +j analysis. For both figures, the lower plot shows the ratio of observed data 
to the combined SM expectation, with the inner band representing the range of statistical uncertainty and 
the outer band (visible only in the highest Rj bin in the right-hand figure) the combined statistical and 
systematic uncertainty. Events outside the range of the displayed region are included in the highest-valued 
bin. 
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A summary of the selection requirements specific to each of the diphoton SRs is presented in Table 2, 
to SRj' b , SRjif SRj ' 1 and SR^ in Table 3, and to the two photon+f SRs in Table 4. After all selection 
requirements, the numbers of events remaining in the various signal regions are 0 (SR^ ( , SR" h ) , 5 
(SR" L ), 1 (SR"_ h ), 12 (SR^ b ), 2 (SR*\ SR[ j , SR^ j ), 16 (SR*) and 10 (SR*). 


Table 2: Enumeration of the requirements defining the four SRs developed for the diphoton signature 
search. 


Signal Region 

SRf-L 

sr£ h 

SR w-l 

ocTT 

^ k w-h 

Number of photons ( Ej [GeV]) 

> 1 (> 75) 

> 1 (> 75) 

> 1 (> 75) 

> 1 (> 75) 

£.™ ss [GeV] 

> 150 

> 250 

> 150 

> 200 

Hj [GeV] 



> 600 

>400 

m e ff [GeV] 

> 1800 

> 1500 



A0 m in(jet, /i' | n ' ss ) (Number of leading jets) 

> 0.5 (2) 

> 0.5 (2) 

> 0.5 (2) 

> 0.5 (2) 

A0 min (y,£7 ss ) 


>0.5 


>0.5 
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Table 3: Enumeration of the requirements defining the four SRs developed for the photon+Z? and photon +j 
signature searches. 


Signal Region 

SRp 

SR? 

sr l 

SRi 

Number of photons (Ej [GeV]) 

> 0 (> 125) 

> 0(> 150) 

1 (> 125) 

1 (> 300) 

E™ iss [GeV] 

> 100 

>200 

>200 

> 300 

Hj [GeV] 


> 1000 


> 800 

Number of jets (number of b-jcts) 

2-4 (> 1) 

> 3 (> 0) 

> 3 a 

> l a 

Number of leptons 

0 


0 

0 

M bh [GeV] 

75 - 150 




v i?miss 

Mj T [GeV] 

> 90 

>90 



A 0 m in(jet, E.|" ss ) (number of leading jets) 

> 0.3 (2) 

> 0.3 (4) 

>0.4 (2) 

> 0.4 (2) 

R* 



<0.85 


^ 0 min(j^’ y) 




< 2.0 


a For S R j’ 1 and S R '- 1 . the two leading jets are required to have p-y > 100 and p-y > 40 GeV, respectively. 


Table 4: Enumeration of the requirements defining the two SRs developed for the photon+f signature 
search. 


Signal Region 

SR y f 

SR? 

Number of photons (Ej [GeV]) 

> 0 (> 125) 

> 0(> 125) 

E™ ss [GeV] 

> 120 

> 120 

//f s [GeV] 

< 100 

< 100 

Number of leptons 

> 0 (e) 

>0 (p) 

1 M ey - M z \ [GeV] 

(> 15) 


p zpmiss 

M t t [GeV] 

> 120 

> 120 

A R{£,y) 

>0.7 

>0.7 
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7 Background estimation 


Backgrounds to the various SRs arise from a number of sources, including processes such as radiative 
vector boson production that generate real photons in combination with energetic neutrinos, as well as 
events in which one or more energetic jets or electrons are misidentified as a photon. While these sources 
contribute generically to all four signatures explored in this study, the differing definitions of each of the 
associated SRs lead to, in many cases, significant differences in the manner in which the contributions 
of these various background sources are estimated. In the following, the methodology of the background 
estimation for each of the four experimental signatures is discussed, and the resulting background estim¬ 
ates, broken down by source, are tabulated. For the estimation of background contributions that rely upon 
MC simulation, either directly or through the estimation of “transfer factors” relating the background con¬ 
tent of control regions to that of corresponding SRs, the effect of MC modeling uncertainties have been 
considered; in general, these uncertainties are found not to be dominant contributions to the overall uncer¬ 
tainty in the background estimates. Background models are confirmed in validation regions (VRs) with 
selection criteria closely related to those of the corresponding SR, but with one or more selection criteria 
modified to suppress the contribution of possible GGM signal to the VR. 

7.1 Backgrounds to the diphoton analysis 

Backgrounds from SM contributions to the four diphoton SRs are grouped into three primary components. 
The first of these, referred to as “QCD background,” arises from a mixture of processes that include yy 
production as well as y + jet and multijet events with at least one jet misreconstructed as a photon. The 
second background component, referred to as “EW background,” is due to W + X (here “ X ” can be any 
number of jets, accompanied by no more than one photon; the two-photon case is treated separately) and 
tt events, with a smaller contribution arising from Z + X events. These events tend to include final-state 
neutrinos that produce significant £™ 1SS . In both cases, EW background events entering the signal regions 
generally have at least one electron misreconstructed as a photon. The QCD and EW backgrounds are 
estimated through the use of dedicated control samples of data events. 

The third background component, referred to as “irreducible,” consists of W and Z bosons produced 
in association with two real photons, with a subsequent decay into one or more neutrinos. For this back¬ 
ground, the W(—> tv) + yy component dominates, and requires corrections to its LO contribution that are 
both large and rapidly varying across the phase space of the W(—> tv) + yy (plus possible additional jets) 
process [73]. Thus a data-driven approach was developed to constrain the W(-» tv) + yy contribution to 
the four SRs. The Z(—> vv) + yy contribution is estimated directly from the MC simulation. 

The QCD background to SR^7 L , SR^ 7 H , SR^ L and SR^i H is expected to arise from events with two 
real, isolated photons (diphoton QCD events) unaccompanied by any additional electroweak bosons, and 
from events with a single real, isolated photon and a jet whose fragmentation fluctuates in such a manner as 
to cause it to be misidentified as a second isolated photon (“photon+jet” QCD events). A contribution from 
dijet QCD events is found to be small and largely incorporated into the photon+jet background estimate. 
To estimate the photon+jet contribution a “QCD control sample” is identified within the diphoton-trigger 
data sample by selecting events for which one photon candidate satisfies the tight selection criterion, while 
the other satisfies the loose but not the tight photon criterion. QCD control sample events containing 
electrons are vetoed to reduce contamination from W —» ev decays. Studies with MC simulated samples 
as well as E™ ss and Hj sideband data suggest that the E | niss distribution of this control sample adequately 
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reproduces the £™ ss distribution of the QCD background in the high-£! [ 11ISS region used for the signal 
selection. A diphoton MC sample is used for the estimation of the diphoton contribution to the QCD 
background. 

The Hj, A0min(jet, C| l " ss ) and A</> m j n (y, £'' | T11SS ) requirements associated with each of the four SRs 
are applied to the QCD control and diphoton MC samples, and the resulting samples are scaled so that 
the combination of the two samples exactly reproduces the number of observed diphoton events (for the 
given SR) in the region 0 < £ , ' | 11ISS < 60 GeV, and with the diphoton MC sample providing a specified 
fraction of the total event count in this region. As suggested by the independent ATLAS H —> yy [74] and 
isolated photon pair cross-section [75] analyses, this fraction is set to 75%, although in this analysis a range 
between 50% and 100% is adopted to reflect the degree of uncertainty in this fraction. The resulting QCD- 
background estimate, for each of the four binolike SRs, is then obtained by summing the scaled number of 
combined QCD control and diphoton MC events with /r' | mss above the minimum requirement for the given 
SR. Additional sources of systematic uncertainty on the QCD-background estimate include its dependence 
on the low-£™ ss region used to scale the diphoton MC and QCD control samples, and the effect of possible 
mismodeling of the A0 m j n (jet, £' | mss ) and A^ m i n (y, E'" ,ss ) distributions of the QCD background by the 
QCD control sample. Including both systematic uncertainty and the statistical uncertainty associated with 
the limited number of events in the QCD control and diphoton MC samples, the result for the QCD 
background and its overall uncertainty is shown in Table 5. 

The QCD-background model is validated by comparing the observed numbers of events to the total 
expected SM background in bins of 300 GeV in Hj for the sideband region 100 < C | n ' ss <150 GeV, for 
which event rates are expected to be dominated by the QCD background. The observed event rate tends to 
be somewhat lower than that predicted by the overall background model, although it is within 1 standard 
deviation of the overall background model uncertainty for all Hj bins. 

The EW background, arising predominantly from W + X and tt events, is estimated via an “electron- 
photon” control sample composed of events with at least one isolated tight photon and one isolated elec¬ 
tron, each with Ej > 75 GeV, and scaled by the probability for such an electron to be misreconstructed 
as a tight photon, as estimated from a “tag-and-probe” study of the Z boson in the ee and ey sample. The 
electron-to-photon scale factor varies between 1.9% (0 < I77I < 0.6) and 3.7% (1.52 < \rj\ < 1.81), since 
it depends on the amount of material in front of the calorimeter. Events with two or more tight photons 
are vetoed from the control sample to preserve its orthogonality to the signal sample. In the case of more 
than one electron, the one with the highest pj is used. Including systematic uncertainties of ±25% each, 
associated with a possible pj dependence of the scale factor and a possible overlap between the QCD and 
EW background estimates, leads to the estimates for the EW background to the four diphoton SRs shown 
in Table 5. 

The irreducible background is composed of two distinct components: diphoton production in asso¬ 
ciation with either aWorZ boson. The latter contribution is relatively small and is sufficiently well 
understood to allow the use of the MC simulation, with a total cross section scaled to that of Ref. [58], 
to directly estimate the Z(—» vv) + yy contribution to the four SRs. The value of this estimate is shown 
is Table 5; the uncertainty is dominated by a ±50% uncertainty on the Z(—> vv) + yy cross section of 
Ref. [58] that arises from the variation of the factorization and renormalization scales used to quantify the 
uncertainty due to missing higher-order processes. 

The W (—» tv) + yy background to the four SRs is estimated using a lepton-diphoton (tyy) CR. To 
enhance the contribution of W(—> tv) + yy and ensure that the tyy CR is exclusive of the four SRs, the 
photon pj requirement is lowered to 50 GeV and a requirement of 50 < E™ ss <150 GeV is imposed. 
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To ensure that the CR sample arises from the same region of the W(—> tv) + yy process phase space 
as the expected background, a further requirement that the transverse momentum of the tyy system be 
greater than 100 GeV is imposed. A total of seven events is observed in the CR, for which MC simulation 
suggests that 2.2 arc expected to arise from SM sources other than \\'(—> tv) + yy. When setting limits 
on contributions from new physics in the four SRs, a simultaneous fit to the CR and the signal region 
under study is performed, allowing both the signal and W(—> tv) + yy contributions to float to their best-fit 
values. When setting model-dependent limits, the fit also takes into account a possible signal contribution 
to the tyy CR, which can be significant for the electroweak-production models in the case that thejfj mass 
is light. In the limit that no GGM signal contributes to the tyy control region, an enhancement factor of 
2.3 must be applied to the W (—» tv) + yy MC sample to achieve agreement between the MC simulation 
and data in the tyy control region. The resulting W(—> tv) + yy-background estimate in each of the 
four SRs, under the assumption that there is no signal contribution to the tyy CR, is shown in Table 5; 
the uncertainty is dominated by that of the limited number of events in the tyy CR. Also shown is the 
combined background estimate, including uncertainty, from all four sources. 

7.2 Backgrounds to the photon+/> analysis 

For both SR' b and SR^ b , which include a requirement of at least one /;-jet, backgrounds arise from two 
predominant sources: from leptonic decays of real or virtual W bosons accompanied by the production of 
ft-quark pairs, including those arising in it events [“W(—» Tv)” backgrounds]; and from events containing 
no electroweak bosons or top quarks (QCD backgrounds). VT(—> IV) background events are further classi¬ 
fied according to the origin of the high-energy isolated photon. Contributions from W(—» tv) backgrounds 
for which the photon arises from the misidentification of an electron are estimated via a control sample 
for which the photon requirement is replaced by an electron requirement, scaled by an electron-to-photon 
misidentification probability; this approach is similar in nature to that of the diphoton analysis. Estimates 
of this component of the background to SRj‘ b and SR^ b are shown in Table 6 ; the quoted uncertainty arises 
from the limited number of events in the control sample, as well as systematic uncertainty associated with 
the possible p-\ dependence of the electron-to-photon misidentification-rate scale factor. 

Contributions from W (—» tv) backgrounds for which the photon is real, or for which the photon arises 
from a misidentified jet or r lepton, are estimated via lepton-enriched CRs that constrain the normalization 

Table 5: The expected and observed numbers of events for the four diphoton signal regions. The quoted 
errors are the combined statistical and systematic uncertainties. 


Signal Regions 

sr£l 

sr£h 

SR w-l 

c R ry 

jr w-h 

Expected background events 
QCD 

EW 

0.06$g 
o oo +0 - 24 

u - uu -o.oo 

0.02 ± 0.02 

0-06!j$ 

o oo +0 - 24 
u - uu -o.oo 

0.0 ± 0.0 

2 04 +0 - 82 
-Q.75 

0 32+ 0 - 45 
u " 5 -0.32 

0.64 ± 0.27 

1 oi +0 - 48 

1 u1 -0.42 

0 22 +033 
-0.22 

0.13 + 0.08 

{W — > tv)yy 

0.04 ± 0.02 

0.05 ± 0.04 

1.01 ±0.62 

0.53 ± 0.34 

(Z — > vv)yy 

0.00 ± 0.00 

0 . 01 + 0.01 

0.07 ± 0.04 

0.13 + 0.07 

Observed events 

0 

0 

5 

1 
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of MC samples used to simulate contributions from these two sources. Separate control regions CRj Lp and 
CR^ P are defined for the low- and high-neutralino-mass SRs by requiring a lepton in addition to the 
requirements already imposed to define the SRs. In addition, in order to increase the number of events 
in the CR, the Z?™ ss requirement is reduced, the M/,/, requirement is removed (for the SRj' b analysis), 
and the H\ requirement is relaxed (for the SR^ b analysis). Events in these two CRs are expected to be 
dominated by tt, tty and Wy production, as is expected for the corresponding background contributions 
to the SRs, and any overlapping phase space is subtracted as part of the background estimation. Including 
all SM sources, a total of 14.5 (58.0) events arc expected in the CR^ P ( CR[*| p ) control regions, to be 
compared to an observation of 18 (61) events. Scaling the combined SM MC samples by these ratios of 
data to expectation, after having subtracted the contributions estimated by other techniques, yields the SR 
background estimates shown in Table 6. It is found that the data-to-expectation scale factor is somewhat 
dependent upon the requirements used to define the lepton-enriched CRs; these variations are included in 
the systematic error on the resulting SR background prediction. 

The QCD background is estimated via the definition of a two-dimensional signal- and control-sample 
grid (the “ABCD” method). For the SRj b analysis, three control samples are defined by requiring only 
a single tagged 6-jet, by requiring that £™ ss <75 GeV, or by requiring both of these SR modifications. 
For the SR'' b analysis, three similar control samples are defined by requiring that no jet be identified as 
a 6-jet, by requiring that i?™ ss <150 GeV, or by requiring both of these SR modifications. A transfer 
factor is calculated by taking the ratio of the number of events with only the C | T " SS requirement changed 
to the number of events with both the E™ ss and 6-jet requirements changed. Assuming that the relaxation 
of the 6-jet requirement is uncorrelated with the relaxation of the Ej lss requirement, the number of QCD- 
background events in the SR can then be estimated by scaling, by this transfer factor, the number of 
events with only the 6-tag requirement changed. This scaling is done only after subtracting the number of 
events expected to come from sources other than those that produce QCD-background events from each 
of the control samples. To avoid the biasing effects of possible correlations between the relaxation of the 
6 -jet requirement and the /7. l [ niss requirement, for the SR' b (SR^ b ) analysis events are binned in and 
weighted bin by bin in the ratio of the number of events in the 2-tag (1-tag) region to the number of events 
in the 1-tag (0-tag) region in the y+jct MC sample. The resulting estimate of the small expected QCD 
background in the two SRs is shown in Table 6, with the systematic uncertainty dominated by the limited 

Table 6: The expected and observed numbers of events for the two photon+6 signal regions. The quoted 
errors are the combined statistical and systematic uncertainties. 


Signal Regions 

SR? 

SR? 

Expected background events 

18.8 ±5.3 

3.82 ± 1.25 

e —> y 

3.2 ± 0.4 

0.18 ±0.08 

W(-» tv) 

12.6 + 4.9 

3.35 ± 1.05 

QCD 

2.3 ±2.1 

0.00 ± 0.65 

Z — > vv 

0.8 ± 0.4 

0.29 ±0.15 

Observed events 

12 

2 
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number of events to which the scale factor is applied. 

An additional background due to the production of a Z boson that decays into two neutrinos, in associ¬ 
ation with a photon and a 6-jet, is estimated directly from the MC simulation, and is tabulated in Table 6. 
For this final contribution, a 50% scale error is assumed for the overall rate of production for this process. 
The combined background from all expected sources is also shown in Table 6. 

For both photon+6 SRs, the background model is validated in four VRs, defined for SR| h by requiring 

2jmiss 

75 < £™ ss <100 GeV, by reversing the requirement, by requiring Mj,’ T <90 GeV, or by requiring 
ArAminijct, E' | mss ) < 0.3, respectively. Since no M^, requirement is made for SR^ b , the second validation 
region is instead defined by changing the Hj requirement to 500 < Hj < 1000 GeV. The observed 
numbers of events in the VRs arc consistent with the predictions of the overall background model. 

7.3 Backgrounds to the photon-i- j analysis 

Backgrounds to the photon+j analysis are expected to arise both from events with real photons as well 
as events for which an electron or a jet is misidentified as a photon. The former source is expected to 
receive contributions from events for which a W/Z boson, a single top quark, or a tt pair is produced in 
association with a real photon, with neutrinos in the subsequent weak decays of these produced states 
providing significant £™ ss (Wy, Zy and tty background). Events with real photons can also contribute to 
the background to the photon +j analysis when significant Z?™ ss arises from instrumental sources (QCD 
background). The Wy, tty and QCD backgrounds are estimated by scaling a corresponding MC sample 
to match the observed event count in a corresponding CR enriched in the given background process but 
otherwise kinematically similar to the corresponding SR. The MC simulation is then used to provide an 
estimate of the expected background in the SR‘ >J and SR^J SRs. Smaller contributions from single-top+y 
and Zy are estimated directly from the MC simulation. 

The QCD-background CR is defined by changing the SR' 1 and S R j^ 1 E™ ,ss requirements to instead 
select events with £™ ss < 50 GeV, but leaving all other selection requirements unchanged, providing 
a region dominated by real photons arising from radiative QCD processes. The VV'y-backg round CR is 
defined by requiring, in addition to the other SRj' J and SR^ requirements, that there be a single identified 
isolated lepton (electron or muon) and no 6-jet in the event. The tty-background CR is defined similarly, 
but requires instead at least one 6-jet. In both cases, in order to increase the number of events in the CR 
the £' | '" ss requirement is changed to 100 < E™ 1SS < 200 GeV. The event counts in the resulting QCD, Wy 
and tty CRs are used to scale the y+jet, Wy and tty MC samples, respectively, after applying a selection 
identical to that of the corresponding CR. The scale factors are determined in a simultaneous fit to all 
CRs, taking into account mutual cross contamination between the different backgrounds. Estimates for 
the contributions of all three of the real-photon backgrounds are shown in Table 7. Systematic uncertainty 
on the scale factor is dominated by the theoretical uncertainties on the relevant MC samples, related in 
turn to the PDF choice and the renormalization and factorization scales. 

As in the other analyses, backgrounds from events for which electrons are misidentified as photons are 
estimated by identifying a control sample of events through the application of a set of selection require¬ 
ments that are identical to those of the given SR, but with a requirement that the event have an electron that 
replaces the required photon. The estimate of the background in the SR (SR^ J or SR^jj) is then, as in the 
other analyses, derived by scaling each event in the control sample by an //-dependent electron-to-photon 
misidentification factor. The resulting background estimates are displayed in Table 7. 
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Finally, the contribution of a background due to events for which the selected photon arises from 
the misidentification of a jet is estimated by determining the jet-to-photon misidentification rate from 
the observed isolation-energy distribution of energy in a cone of size A R = 0.2 surrounding the energy 
deposition in the calorimeter associated with the photon. The isolation-energy distribution for real photons 
is modeled with electrons from Z boson decays, while that of misidentified jets is modeled with a sample 
of events for which there is a “pseudophoton.” A pseudophoton is defined to be an object that passed all 
loose photon selection requirements, as well as all tight photon selection requirements except one or more 
from a set of four that relate to the shape of the deposition in the finely granulated front portion of the EM 
calorimeter. The fraction of misidentified jets within the tight, isolated photon sample is determined with 
a control sample composed of events with tight, isolated photons with pj >125 GeV, as well as a relaxed 
E'| llss requirement of 50 < T , ' | 11ISS < 150 GeV and an intermediate requirement of Hj > 600 GeV. The 
photon isolation-energy distribution of this control sample is fit to establish the relative amounts of these 
two sources (real photons and misidentified jets), with the misidentification fraction taken to be the relative 
integrals of the isolation-energy distributions of the misidentified and total contributions in the region for 
which the isolation energy is less than 5 GeV. The estimation of the jet-misidentification background 
in each signal region and control sample (as well as for the validation regions described below) is then 
obtained by scaling the observed number of events in each region or sample by the jet-misidentification 
factor. The number of misidentified jets is then parametrized as a function of F™ ss by fitting the E™ ss 
dependence of the estimated misidentified-jet contribution in the range £™ ss < 200 GeV. The estimates 
in the S R j 1 and SRj'j 1 signal regions are then extracted by integrating the fit function over the relevant 
E™ ss range. The result for the contribution of the jet-misidentification backgrounds for each SR is shown 
in Table 7. Systematic uncertainties arise due to the uncertainties in the combined fit used to derive the 
misidentification factor, for which the parameters of the signal and background templates are allowed to 
vary within then - uncertainties, and from the uncertainties of the extrapolation fit used for the estimation 
of the SR contamination. 

The background model is validated by comparing expected and observed event rates in several VRs. 
For SR”, this includes three VRs for which the A0 m i n (jet, £™ ss ) is reversed, £™ ss is required to be within 
an intermediate range of 75 < £™ 1SS <150 GeV, and for which the R 4 y requirement is reversed. For SR^ J 
two VRs are made use of, including one for which A0 m j n (jet, E'™ ss ) is reversed and another that requires 
that 400 < Hj < 800 GeV. Good agreement is observed between the number of expected and observed 
events in all five VRs. 

7.4 Backgrounds to the photon+7 analysis 

Backgrounds to the photon+7 analysis (SR^ and SRj/j are expected to arise primarily from events with 
hard photons produced in association with electroweak bosons (Wy or Zy) and top quarks {tty), and 
events containing W bosons or semileptonically decaying top quarks for which an accompanying jet is 
misidentified as a photon (jet-to-photon events). Fesser contributions are expected to arise from tt events 
and events containing two electroweak bosons that produce two final-state leptons, one of which is an 
electron that is subsequently misidentified as a photon. As in the other analyses, data-driven techniques 
making use of CRs similar to but exclusive of the SRs, or control samples appropriate for assessing jet- 
to-photon and electron-to-photon misidentification rates, are used to estimate or constrain the primary 
backgrounds, while lesser backgrounds are estimated directly from MC simulation. 

The most prevalent background in the photon+7 sample is expected to arise from Wy events. A Wy 
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CR is defined by requiring an isolated electron or muon, and by requiring in addition that 45 < £™ ss < 

£ ^pmiss p 

100 GeV and 35 < M T T < 90 GeV, but otherwise requiring that the sample satisfy the SR;/ and 
SRj’f criteria. Transfer factors relating the number of events observed in the Wy CR to the number of 
Wy events expected in the SRs are estimated, separately for the electron and muon contributions, from 
the Wy MC simulation. Systematic uncertainties on the resulting VV'y-background estimate for the two 
SRs arise from the scale and PDF uncertainties associated with the transfer factors. A somewhat lesser 
contribution from tty events is estimated directly from the MC simulation, with uncertainties arising from 
imprecise knowledge of the strong-interaction scale and the rate of final-state photon production into the 
acceptance of the SRs. A smaller background contribution from Zy events is estimated directly from the 
MC simulation, with an uncertainty of ±50% assumed for the production rate into the region of the Zy 
phase space that populates the photon+f SRs. 

As in the other analyses, a potentially sizable contribution to the photon ±^ sample arises from jet- 
to-photon misidentification. The contribution of SR events arising from jet-to-photon misidentification is 
estimated by exploring the isolation-energy distribution of events in an extended Wy control sample for 
which the requirement on isolation energy has been removed. Isolation-energy distribution templates for 
true photons and for jets misidentified as photons are developed in the manner described for the photon +j 
analysis. A tit is then performed on the isolation-energy distribution of the extended Wy control sample 
to estimate the number of events in the isolated (isolation energy less than 5 GeV) Wy CR that arise from 
jets misidentified as photons. A scale factor is defined as the ratio of the estimated number of events in 
the isolated Wy CR arising from misidentified jets to that expected from the VF+jcts and semileptonic 
it MC simulations. A data-driven estimate of the number of events arising from misidentified jets in 
SR;^ and S Rj/ is then derived by multiplying the number of such events expected from the combination 
of the VV'+jcts and semileptonic ft" MC simulations by this scale factor. Because the MC simulation is 
relied upon to propagate the background estimate from the control sample into the SR, uncertainties on 
the jet-to-photon misidentification background arise due to imprecise knowledge of the proton PDFs and 
strong-interaction scale. An additional uncertainty is assigned based on the difference between the scale 
factors determined for the separate electron and muon samples. 

A final significant source of background is expected to arise from it events, single-top events, and 
events containing two electroweak bosons that produce two final-state leptons, one of which is an electron 
that is subsequently misidentified as a photon. The contribution from these backgrounds is estimated 
from MC simulation, applying a correction based on the relative electron-to-photon misidentification rate 
between data and MC simulation. This correction is determined from Z —> e + e~ events as described 
above for other analyses. In addition to the uncertainty in the measurement of the misidentification rate, 
uncertainties in the estimate arise from PDF and scale uncertainty in the tt production process as well as 
an assumption of a ±50% uncertainty in the rate of single-top and diboson production. 

All other sources of background, including those from Z+jet, y+jet and yy production, are expected 
to contribute only minimally to the total SR backgrounds. In particular, a potential background from 
y±jet events arising from jet-to-lepton misidentification is estimated using a matrix method (as described 
in Ref. [76]) making use of a control region incorporating nonisolated lepton candidates, and is found 
to contribute 0.1 events to the overall background estimate for each of the SRj f and SR^ samples. A 
summary of the resulting background estimates for the SRj f and SRjf SRs is shown in Table 8, broken 
down by source. 

The background model is validated for each SR by comparing expected and observed event rates in 
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two VRs. An M\ VR is defined by relaxing the M T 1 requirement to 35 < M T T <90 GeV; to 
increase the number of events in this VR, the £' | mss requirement is also relaxed to £™ ss >100 GeV. A 

£ rmiss 

Ej lss VR is defined by relaxing the fs™ ss requirement to 45 < £'™ 1SS <100 GeV while leaving the M T T 
requirement unchanged. For the electron and muon channels combined, the number of events in the E' | 11ISS 
VR is observed to be somewhat less than that expected for the background model, although still within 2 
standard deviations of the combined statistical and systematic uncertainty. Good agreement is found for 
the Mj VR. 

8 Signal efficiency and systematic uncertainty 

GGM signal acceptances and efficiencies are estimated using MC simulation for each simulated point 
in the gluino-bino, wino-bino, higgsino-bino, and wino-NLSP parameter spaces, and vary widely across 
the regions of these spaces relevant to establishing the limit contours presented below. The product of 
acceptance times efficiency tends to be greatest (10%-25%) when the masses of both the produced and the 
NLSP states are largest, leading to large amounts of both visible energy and missing transverse momentum 
that would clearly distinguish signal from background events. However, for the lower accessible mass 
scales associated with electroweak production, and particularly for the case of a low-mass NLSP, the 
product of acceptance times efficiency can be significantly smaller. For example, for the region relevant 
to establishing limits at low values of p, the efficiency of the SR” analysis is less than 0.1%, leading to a 
relatively modest lower limit on the mass of produced SUSY states. 

Making use of a bootstrap method [77], the efficiencies of both the single photon and diphoton triggers 
are determined to be greater than 99%, with an uncertainty of less than 1%. 

The reconstruction efficiency for tight, isolated photons is estimated with complementary data-driven 
methods [78]. Photons identified kinematically as having come from radiative decays of a Z boson (Z —> 
t + l~y events) are used to study the photon reconstruction efficiency as a function of p-\ and //. Independent 
measurements making use of a tag-and-probe approach with Z —> ee events, with one of the electrons 
used to probe the calorimeter response to electromagnetic depositions, also provide information about the 
photon reconstruction efficiency. For photons with Ej > 75 GeV, the identification efficiency in the range 
0 < |? 7 l < 1.81 is greater than 95%; for the range 1.81 < \r/\ < 2.37 the efficiency is approximately 90%. 
The uncertainty in the efficiency also varies with \rj\, and lies between ±(l-2)% 

The isolated electron efficiency is also estimated using tag-and-probe methods, making use of samples 
of Z ee and J/iJj —> ee events as described in Refs. [62, 63]. The efficiency and its uncertainty are 
estimated as a function of electron p-\ and //. leading to an overall uncertainty of ±1.0% on the efficiency 
of the photon+f analysis, the only analysis that explicitly requires an electron. The muon identification 
uncertainty, estimated as described in Ref. [65], is found to contribute an uncertainty of only 0.4% on the 
efficiency of the photon+Z’ analysis. 

In portions of the GGM parameter space, uncertainties that vary across the parameter space dominate 
the systematic uncertainty on the signal acceptance times efficiency. These model-dependent uncertainties 
include those due to uncertainties in the photon, electron and jet-energy scales, the fi-jet tagging efficiency, 
and the “pileup” uncertainty arising from the modeling of additional interactions in the same or nearby 
bunch crossings. 

The electron and photon energy scale is determined using samples of Z —> ee and .//i/z —> ee 
events [79], both of whose masses are known precisely and thus provide an accurate calibration signal 
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Table 7: The expected and observed numbers of events for the two photon +j signal regions. The quoted 
errors arc the combined statistical and systematic uncertainties. 


Signal Regions 

SR l 
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Expected background events 
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Table 8 : The expected and observed numbers of events for the two photon+f signal regions. The quoted 
errors arc the combined statistical and systematic uncertainties. The contribution from the Zy process 
arises from events for which one of the leptons from the Z —* ( + (~ decay is either missed or badly 
mismeasured. The likelihood of this occurring is significantly greater for muons than electrons. 

Signal Regions 

SR l ( 

SR f 

Expected background events 

10.5 ± 1.4 

14.1 ± 1.5 

Wy 

6.7 ± 1.2 

8.8 + 1.3 

tty 

1.4 ±0.6 

1.7 ±0.7 

Zy 

0.0 ± 0.0 

1.2 ± 0.6 

Jet —> y 

1.5 ± 1.0 

1.2 ±0.7 
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0.7 ±0.2 

0.8 ±0.3 

Other sources 

0.3 ±0.1 

0.4 ± 0.2 

Observed events 

16 
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for determination of the electromagnetic calorimeter response. Uncertainties arise from imprecise know¬ 
ledge of the material burden between the IP and the face of the EM calorimeter. The muon energy scale 
and uncertainty are similarly estimated with calibration samples of Z —> pp, T —> pp and J/tp —> pp 
events [65]. 

The jet-energy scale is established via the propagation of single-particle test-beam measurements of 
the calorimeter response through simulations of jets arising from pp collisions [67, 80]. The jet-energy 
scale uncertainty is constrained by the study of momentum imbalance in dijet events [81], as well as from 
an assessment of the effect of uncertainties in the modeling of jet properties with MC simulations, and 
from uncertainties in the modeling of the varying response to differing jet flavor composition. 

Uncertainties in the values of whole-event observables, such as E™ ss and Hj, arise from uncertainties 
on the energy of the underlying objects from which they are constructed. In addition, the £™ 1SS observable 
receives a contribution from calorimetric energy deposits not associated with any of the reconstructed 
objects in the event. Uncertainties on the energy scale of these unassigned contributions are found to 
contribute negligibly to the overall uncertainty on the value of the ZT | mss observable. 

The uncertainty due to pileup is estimated by varying the distribution of the number of interactions 
per bunch crossing overlaid in the simulation by ±10%. The uncertainty on the ^-tagging efficiency in the 
MC simulation is estimated from measurements of dedicated heavy-flavor calibration data samples. 

In the regions of GGM parameter space relevant for establishing the exclusion limits discussed be¬ 
low, and excepting MC statistical uncertainty, the quadrature sum of the individual sources of systematic 
uncertainty on the signal reconstruction efficiency for the diphoton, photon+b and photonanalyses is 
of order 10%. For the photon +j analysis the systematic uncertainty is somewhat larger—approximately 
20 %— due to an increased sensitivity to the jet-energy scale and resolution associated with the multiple-jet 
requirement. 

9 Results 

An accounting of events observed in each SR is shown in Table 9, along with the size of the expected SM 
background. Comparisons of the C [ niss distribution between signal and expected background is shown for 
several different SRs in Figs. 5-7. No evidence for physics beyond the SM is observed in any of the SRs. 
The largest excess relative to the expected background is observed for the SR^J_ L analysis; considering 
both statistical and systematic uncertainty, and assuming that all observed events are from SM sources, 
an observation of five or more events over an expected background of 2.04^ ^ represents an upward 
fluctuation with a probability of occurrence of approximately 6%. 

Based on the numbers of observed events in the ten SRs and the background expectation shown in 
Table 9, 95% confidence-level (CF) upper limits are set for each SR on the number of events from any 
scenario of physics beyond the SM, using the profile likelihood and CL S prescriptions [82]. Uncertainties 
on the background expectations are treated as Gaussian-distributed nuisance parameters in the maximum- 
likelihood fit. Assuming that no events due to physical processes beyond those of the SM populate the 
various CRs used to estimate SR backgrounds, observed 95% CF limits on the number of such events vary 
between 3.0 (for the SR^7 L and SR^7 H SRs) and 14.2 (for the SRj f SR). Taking into account the integrated 
luminosity of 20.3±0.6 fb _1 these number-of-event limits translate into 95% CF upper limits on the visible 
cross section for new physics, defined by the product of cross section, branching fraction, acceptance and 
efficiency for the different SR definitions. Correspondingly, the observed visible cross-section limits vary 
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Table 9: Summary of the number of events expected from SM sources (N^), and the observed number 
of events (AobsX for each of the ten SRs. Also shown is the derived model-independent 95% CL limit 
(S*h on the number of possible events from new physics, as well as both the observed ((ecr) 9 ^) and 
expected ((£cr)exp) 95% CL limit on the visible cross section from new physics. Due to the discrete nature 
of the number-of-observed-events likelihood distribution in background-only pseudoexperiments, when 
both the expected number of background events and its uncertainty arc close to zero the expected limit is 
dominated by the case of zero observed events. This leads to a very narrow one-standard-deviation range 
for the expected limit for SR^7 L and SR^7 H . 
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Figure 5: Distribution of the missing transverse momentum E™ ss for the sample surviving all requirements 
of the SR^ H (left) and SR^_ L (right) selection except the E'™^ requirement itself. Overlain are the 
expected SM backgrounds as a function of £ , ' [ 111SS , separated into the various contributing sources. Also 
shown are the signal expectations for the (in^,. m~») = (600,100) GeV and (mp, m~<>) = (600,500) GeV 
models. The lower plots show the ratio of observed data to the combined SM expectation. For these plots, 
the inner band represents the range of statistical uncertainty while the outer band represents the combined 
statistical and systematic uncertainty. Events outside the range of the displayed region are included in the 
highest-valued bin. 
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Figure 6: Z7 | niss distribution for the sample surviving all requirements of the SRj >b (left) and SR^ b (right) 
selection except the /:! | n ' ss requirement itself. Overlain arc the expected SM backgrounds as a function 
of Z?™ ss , separated into the various contributing sources. Also shown are the signal expectations for the 
(in,g,m~o) = (1200,150), (1200,450), and (1200,850) GeV models. The lower plots show the ratio of 
observed data to the combined SM expectation. For these plots, the inner band represents the range of 
statistical uncertainty while the outer band represents the combined statistical and systematic uncertainty. 
Events outside the range of the displayed region arc included in the highest-valued bin. 


27 
























































50 100 150 200 250 300 

£-miss 


-T 


[GeV] 


Figure 7: (Left) i?™ ss distribution for the sample surviving all SR” requirements except the £' | riISS require¬ 
ment itself. Overlain are the expected SM backgrounds as a function of £™ ss , separated into the various 
contributing sources. Also shown are the signal expectations for the two points in the M^~n parameter 
space of the GGM model relevant to the photon-i- j analysis. (Right) £™ 1SS distribution for the combined 
sample of events surviving all SRj f and SR jf requirements except the £. I | 11ISS requirement itself. Overlain 
are the expected SM backgrounds as a function of E ™ ss , separated into the various contributing sources. 
Also shown is the signal expectation for the m^ v = 300 GeV GGM model relevant to the photon+7 ana¬ 
lysis. For both figures, the lower plot shows the ratio of observed data to the combined SM expectation, 
with the inner band representing the range of statistical uncertainty and the outer band the combined stat¬ 
istical and systematic uncertainty. Events outside the range of the displayed region are included in the 
highest-valued bin. 
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between 0.15 and 0.70 fb. 

By considering, in addition, the value and uncertainty of the acceptance times efficiency of the se¬ 
lection requirements associated with the various SRs, as well as the NLO (+NLL) GGM cross sec¬ 
tions [33-37], which vary steeply with gluino and gaugino mass, 95% CL lower limits may be set on 
the masses of these states in the context of the various GGM scenarios explored in this study. For the 
diphoton, photon +6 and photon-i -j analysis, the SR with the best expected sensitivity at each simulated 
point in the parameter space of the corresponding GGM model(s) is used to determine the degree of ex¬ 
clusion of that model point. For the photon+f analysis, the 95% CL exclusion limits are derived from the 
combined likelihood of the electron and muon channels, taking into account the correlation between the 
systematic uncertainty estimates in the two channels. 

For the diphoton analysis, SR^7 r is expected to provide the greatest sensitivity to the gluino-bino 
model for bino masses above 800 GeV and SR^7 L for bino masses below this. For the wino-bino model, 
the similar transition point between the use of SR^J L and SR^ H is found to be at 350 GeV. The resulting 
observed limits on the gluino and wino masses are exhibited, as a function of bino mass, for the diphoton 
analysis gluino and wino production models in Figs. 8 and 9, respectively. For the puipose of establishing 
these model-dependent limits, for all four diphoton SRs both the normalization of the W (—» (v) + yy- 
background estimate and the limit on the possible number of events from new physics arc extracted from 
a simultaneous fit to the SR and W (—> tv) + yy control region, although the signal contamination in the 
W(—> (:v) + yy control sample is appreciable only for the low-bino-mass region of the wino-bino parameter 
space. Also shown for these two figures, as well as for the following two figures (Figs. 10 and 1 1 ), are the 
expected limits, including their statistical and background uncertainty ranges, as well as observed limits 
for SUSY model cross sections ±1 standard deviation of theoretical uncertainty from their central value. 
Conservatively choosing the -1 standard-deviation observed contour, 95% CL lower limits of 1290 GeV 
and 590 GeV are set by the diphoton analysis on the value of the gluino or wino mass, respectively, for 
any value of the NLSP bino mass less than that of the gluino (wino) mass. 

Due to the discrete nature of the number-of-observed-events likelihood distribution in background- 
only pseudoexperiments, when both the expected number of observed events and its uncertainty are close 
to zero the expected limit is dominated by the case of zero observed events. This leads to a very narrow 
one-standard-deviation range for the expected limit, as observed for the expected-limit contour displayed 
in Fig. 8 . In addition, because the observed number of events is very close to the expected number of 
events for SR^ H and SR(/ l , the expected and observed limits are nearly identical in Fig. 8 . 

For the photon +6 analysis, limits are set in the two-dimensional plane of gluino and mass for the 
higgsino-bino GGM model with a negative value of the /i parameter. For NLSP masses near the 95% CL 
exclusion contour, SRj' h is expected to provide greater sensitivity for NLSP masses below approximately 
600 GeV, and so is made use of in this region; above that, SR j' ( h is used to establish the degree of exclusion 
of points in the GGM model space. The resulting observed exclusion contour is shown in Fig. 10 . Again 
choosing the -1 standard-deviation observed contour, in the context of this GGM model a conservative 
lower limit of 1300 GeV is established for the gluino mass over much of the range of the higgsino-bino 
NLSP mass. For NLSP masses above 1000 GeV the sensitivity lessens due to the restriction of the phase 
space for producing an energetic 6 -jet, while for NLSP masses below 600 GeV, the onset of the direct 
production of gaugino states begins to make the analysis insensitive to the value of the gluino mass. 

For the photon+j analysis, limits are set in the two-dimensional plane of the GGM parameters p and 
M 3 for the higgsino-bino GGM model with a positive value of the p parameter. For values of p near 
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GGM: bino-like neutralino, gluino production 



Figure 8: Exclusion limits in the gluino-bino mass plane, using the SR|^ h analysis for m~t> > 800 GeV 
and the SR]7 l analysis for m^o < 800 GeV. Combinations of gluino and bino mass are excluded at 95% 
CL in the area below the unbroken curve. The observed limits are exhibited for the nominal SUSY model 
cross-section expectation, as well as for a SUSY cross section increased and decreased by 1 standard 
deviation of the cross-section systematic uncertainty. Also shown is the expected limit, as well as the ± 1 
and ±2 standard-deviation ranges of the expected limit. 
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Figure 9: Exclusion limits in the wino-bino mass plane, using the H analysis for m^o > 350 GeV 
and the SR^_ L analysis for m~o < 350 GeV. The vertical axis represents wino mass while the horizontal 
axis represents bino mass. The observed limits are exhibited for the nominal SUSY model cross-section 
expectation, as well as for a SUSY cross section increased and decreased by 1 standard deviation of 
the cross-section systematic uncertainty. Also shown is the expected limit, along with its ±1 standard- 
deviation range. The discontinuity at mg o = 350 GeV is due to the switch between the use of the SR^_ L 
and SR^_ H analyses, the former of which exhibits a small excess of observed events relative to the expec¬ 
ted SM background. 
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GGM: higgsino-like neutralino, p.<0 



Figure 10: Exclusion limits in the gluino-neutralino mass plane, for the higgsino-bino GGM model with 
H < 0, using the merged (see text) SR^ b and SR^ b analyses. The observed limits are shown for the nominal 
SUSY model cross-section expectation, as well as for a SUSY cross section increased and decreased by 
1 standard deviation of the cross-section systematic uncertainty. The expected limit is also shown, along 
with its ±1 cr range. For NLSP masses below approximately 450 GeV, the onset of the direct production 
of gaugino states makes the analysis insensitive to the value of the gluino mass. 
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the 95% CL exclusion contour, SR£ j is expected to provide a greater sensitivity for NLSP masses below 
approximately 900 GeV, and so is made use of in this region; above that, SR^ J is used to establish the 
degree of exclusion of GGM model-space points. The resulting observed exclusion contour is shown in 
Fig. 11 . Again choosing the -1 standard-deviation observed contour, in the context of this GGM model a 
conservative lower limit of 1140 GeV is established for the gluino mass parameter M 3 over much of the 
range of the /u parameter. For values of M 3 close to the value of p for which the gluino mass approaches 
that of the higgsino-bino NLSP, the sensitivity of the analysis lessens due to the restriction of phase space 
for producing multiple high-px jets. 

For the photon-M analysis, a limit is set on the wino mass, the single free parameter of the wino-NLSP 
model. Fig. 12 shows the observed limit on the cross section for wino production in this model, as well 
as the corresponding expected limit with ± 1 and ±2 standard-deviation uncertainty bands. Also shown 
is the cross section as a function of wino mass, with its +1 standard-deviation range. In the context of 
this wino-NLSP model, conservatively choosing the -1 standard-deviation cross-section contour leads 
to an exclusion of GGM winos in the the range 124 < M^ <361 GeV; for M^ < 124 GeV the signal 
contamination in the Wy CR becomes too large to permit a reliable estimate of the Wy background. These 
limits are based on the direct production of the wino NLSP in the limit where squark masses are infinite, 
and are independent of gluino mass. 
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Figure 11: Exclusion limits imposed by the photon +j analysis in the two-dimensional plane of the GGM 
parameters M 3 and p. for the higgsino-bino GGM model with p > 0. using the merged (see text) SR' /J and 
SR|| analyses. The observed limits are shown for the nominal SUSY model cross-section expectation, 
as well as for a SUSY cross section increased and decreased by 1 standard deviation of the cross-section 
systematic uncertainty. The expected limit is also shown, along with its ±1 standard-deviation range. 
Values of M 3 below llOOGeV are excluded for most values of the p parameter, although a significant 
region corresponding to the case for which the gluino mass is close to that of the lightest neutralino 
masses remains unexcluded due to the requirements of one or more jets arising from the gluino decay. 
The top and right axes represent the corresponding values of the lightest neutralino and gluino masses, 
respectively. 
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Figure 12: Contour of exclusion in wino production cross section from the photon+f analysis, as a function 
of the wino mass parameter M^ r . The expected limit is shown along with its ±1 and ±2 standard-deviation 
ranges. Also shown is the total cross section for the production of W pairs, again as a function of 
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10 Conclusion 


Making use of 20.3 fb -1 of pp collision data at \fs = 8 TeV recorded by the ATLAS detector at the 
LHC, a search is performed for photonic signatures of new physics associated with significant E™ ss . Four 
experimental signatures arc explored, each involving at least one energetic isolated final-state photon in 
association with significant E.™ ss , and used to search for evidence for several GGM SUSY scenarios. No 
significant excess of events over the SM expectation is observed in any of the searches and so limits arc 
set on possible contributions of new physics. Model-independent limits arc set on the numbers of events 
from new physics and the associated visible cross section. Model-dependent limits are set on the masses 
of SUSY particles or on mass parameters associated with the various GGM scenario models. 

A diphoton signature is used to explore both strongly and weakly produced SUSY states with a decay 
chain proceeding through a binolike NLSP. In the context of these models, lower limits of 1290 GeV and 
590 GeV arc set on the masses of a degenerate octet of gluinos and a degenerate set of winos, respectively, 
for any value of the bino mass less than the mass of these produced states. A photon-plus-h-jet signature 
is used to search for a scenario in which the GGM NLSP is a higgsino-bino admixture with a roughly 
equal branching fraction to photons and to the SM Higgs boson. In the context of this model, a lower 
limit of 1260 GeV is established for the gluino mass over much of the range of the higgsino-bino NLSP 
mass; for NLSP masses below approximately 450 GeV, the onset of the direct production of gaugino 
states makes the analysis insensitive to the value of the gluino mass. A photon-plus-jet signature is used to 
search for an alternative scenario for which the GGM NLSP is a higgsino-bino admixture with a roughly 
equal branching fraction to photons and to the SM Z boson. In the context of this model, a lower limit 
of 1140 GeV is established for the gluino mass parameter M 3 over much of the range of the higgsino 
mass parameter p. Finally, a photon-plus-lepton signature is used to search for a scenario for which the 
GGM NLSP is a degenerate set of three wino states. Based on the possible direct production of these 
states, in the limit of infinite squark mass GGM winos arc excluded in the range 124 < < 361 GeV, 

independent of the gluino mass. 
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